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Abstract 

We study the effect of scheme dependence upon the NLO QCD analysis of 
the world data on polarized DIS. The reliability of an analysis at NLO is demon¬ 
strated by the consistency of our polarized densities with the NLO transformation 
rules relating them to each other. We stress the importance of the chiral JET 
scheme in which all the hard effects are consistently absorbed in the Wilson co¬ 
efficient functions. 


PACS numbers:13.60.Hb; 13.88-|-e; 14.20.Dh; 12.38.-t 


There has been a major effort in the past decade to obtain reliable information 
about the polarized parton densities in the nucleon, and especially to try to determine 
the degree of polarization of the gluon. Aside from its fundamental interest and its 
potential as a testing ground for QCD, such information is vital for the planning of 
experiments at the RHIC collider, due to come into operation in 1999. 


A great improvement in the quality of the data on inclusive deep inelastic scattering 
of leptons on nucleons has been achieved recently ^ and, concomitantly, several 
detailed NLO QCD analyses of the data have been carried out 1 - 0 . of which the 
most comprehensive are in refs. EM- 


It is well known that at NLO and beyond, parton densities become dependent on 
the renormalization ( or factorization ) scheme employed.]] It is perhaps less well known 
that there are significant differences in the polarized case, as a consequence of the axial 
anomaly and of the ambiguity in the handling of 75 in n dimensions. 


In this paper we study the effect of carrying out the data analysis in different 
schemes, in the MS', AB and what is called the JET scheme, whose importance we 
wish to stress. In the latter all hard effects are consistently absorbed into the Wilson 
coefficient functions. The parton densities in each scheme are, by definition, related to 
each other by certain NLO transformation rules. On the other hand, the evolution 
in each scheme is controlled by the splitting functions (or the anomalous dimensions 
in Mellin n-moment space) relevant to that scheme. Thus, in principle, for any two 
schemes labelled 1 and 2 and for any generic polarized density A/ one has symbolically: 


NLO data analysis in Scheme 1 ^ 

A/i , 

( 1 ) 

NLO data analysis in Scheme 2 

A /2 , 

( 2 ) 

NLO transformation rules on A/i ^ 

A /2 

(3) 


and it is a major test of the stability of the analysis and of the consistency of the theory 
that the results for A /2 in and (^) should coincide. In fact, built into the theoretical 
structure is the feature that the two results actually should differ by terms of NNLO 
order. Hence the degree to which the results agree is a measure of the reliability of 
carrying out an analysis at NLO. Moreover, in the set of schemes we study, the non- 

^Of course, physical quantities such as the virtual photon-nucleon asymmetry Ai{x,Q'^) and the 
polarized structure function gi{x,Q'^) are independent of choice of the factorization convention. 
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singlet and gluon densities are the same in all the schemes. That this feature should 
emerge from the data analysis provides a further test of the stability and reliability of 
the analysis. The importance of an analysis of this kind was first stressed in ref. |P]. 


In the unpolarized case the most commonly used schemes are the MS', MS and 
DIS and parton densities in different schemes differ from each other by terms of order 
<^s(Q^) , which goes to zero as increases. 

There are two signihcant differences in the polarized case: 

ij The singlet densities AE(a:, Q^), in two different schemes, will differ by terms of 
order 

as(Q^)AG(x,Q^) , (4) 


which appear to be of order Og. But it is known |T^, |T^ that, as a consequence of the 
axial anomaly, the hrst moment of the polarized gluon density AG{x, Q^) 


dxAG{x,Q ) oc [as{Q )] 


2m-1 


(5) 


and thus grows in such a way with as to compensate for the factor as{Q^) in (|). 
Thus the difference between AE in different schemes is only apparently of order asiQ"^), 
and could be quite large. 

ii) Because of ambiguities in handling the renormalization of operators involving 
75 in n dimensions, the specification MS does not define a unique scheme. Really 


there is a family of MS schemes which, strictly, should carry a sub-label indicating 


how 75 is handled. What is now conventionally called MS is in fact the scheme due 
to Vogelsang and Mertig and van Neervenen |^, in which the first moment of the 
non-singlet densities is conserved, i.e. is independent of Q^, corresponding to the 
conservation of the non-singlet axial-vector Cabibbo currents. 

Although mathematically correct it is a peculiarity of this factorization scheme 
that certain soft contributions are included in the Wilson coefficient functions, rather 
than being absorbed completely into the parton densities. As a consequence, the hrst 
moment of AE is not conserved so that it is difficult to know how to compare the DIS 
results on AE with the results from constituent quark models at low Q^. 

To avoid these idiosyncrasies Ball, Forte and Ridolh [§] introduced what they called 


the AB scheme, which involves a minimal modihcation of the MS scheme, and for which 
the transformation equations are: 


AT,{x,Q^)ab = AJ:{x,Q^)j^+ Nf 
AG{x,Q^)ab = AG{x,Q^)jjg 


271 


y 


AG{y,Q^) 


MS ’ 


( 6 ) 
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or, in Mellin n-moment space, 


AE(n,Q")^fl = AE(n,Q^)jfs + Nf^^AG(n,Q^)jfs, 

AG(n,Q^)AB = AG(n,Q^)jjs- ( 7 ) 

That AE(n = 1)ab is independent of to all orders follows from the Adler- 
Bardeen theorem iQ. In (P) and (0) iV/ denotes the number of flavours. 

The singlet part of the first moment of the structure function gi 

f'dxg\‘''(x,Q^) (8) 

JO 

then depends on AE and AG only in the combination 


aoiQ^) = AE(1, g2)_ = AE(1 )ab - Nj 


271 


-AG(1,Q' 


(9) 


and the unexpectedly small value for the axial charge Oq found by the EMC [^|, which 
triggered the ’’spin crisis in the parton model” [jT^, can be nicely explained as due 
to a cancellation between a reasonably sized AS(1) and the gluon contribution. Of 
importance for such an explanation are both the positive sign and the large value (of 
order 0{1)) for the hrst moment of the polarized gluon density AG{1,Q‘^) at small 
~ 1 — 10 GeV^. Note that what follows from QCD is that |AG(l,g^)| grows with 
g^ (see Eq. (^)) but its value at small is unknown in the theory at present and has 
to be determined from experiment. 


Although the AB scheme corrects the most glaring weakness of the MS scheme, it 


does not consistently put all hard effects into the coefficient functions. As pointed out 
in [^] one can dehne a family of schemes labelled by a parameter a: 

AE\ ^/'AE\ ^/O /AE\ 

where 



Zgcix] a) = Nf[{2x - l)(a - 1) + 2(1 - x)] , (11) 


in all of which (P) holds, but which differ in their expression for the higher moments. 
(The AB scheme corresponds to taking a = 2). 

Amongst these we believe there are compelling reasons to choose what we shall call 
the JET scheme (a = 1), i.e. 


= 2Nf(l - x) . 


( 12 ) 


This is the scheme originally suggested by Carlitz, Collins and Mueller and also 
advocated by Anselmino, Efremov, Leader and Teryaev in refs. |^, ^.[| In it all hard 
Mhere is misprint in Eq. (8.2.6) of |^. The term ) should be [ln{ ) — !]■ 
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effects are absorbed into the coefficient functions. In this scheme the gluon coefficient 
function is exactly the one that would appear in the cross section for 

PP jet(kr) + jet(-kr) + X , (13) 


i.e., the production of two jets with large transverse momentum k^ and —ky, respec¬ 
tively. 

More recently Muller and Teryaev have advanced rigorous and compelling ar¬ 
guments, based upon a generalization of the axial anomaly to bilocal operators, that 
removal of all anomaly effects from the quark densities leads to the JET scheme. Also 
a different argument by Cheng leads to the same conclusion. (Cheng calls the JET 
scheme a chirally invariant (Cl) scheme.) 

The transformation from the MS scheme of Mertig, van Neerven and Vogeslang to 
the JET scheme is given in moment space by 


AE{n,Q‘^)jET — + 2Nf 

AG{n,Q^)jET = . 


Ois{Q‘^ 


2Tm{n + 1) 


AG(n, Q )j^g , 


(14) 


Note that (p!4D implies that the strange sea As is different in the two schemes. Of 
course, (0) and (p!4D become the same for n = 1. 


The NLO Wilson coefficient functions AC'j-^^(x) and polarized splitting functions 
AP^p{x) (or the corresponding anomalous dimensions A 7 fj^^(n)) for the MS and AB 

|I6| and 0, respectively. The NLO coefficient functions 


schemes can be found in refs. 


and anomalous dimensions in the JET scheme are related to those of the MS scheme 
by mEt 

2Nf 


AC<‘'(n)^ET = ACW(n) 


MS ’ 


Ad(^\n)jET = AC, 


b), 

G 


n. 




MS 


n{n -|- 1) ’ 


(15) 


^lGGi^)jET 


^IqG + 

^7S(’^)ms > 
^lGGi'^)ws ~ 


4X/ 

n(n -|- 1) 



1 


4iV/ 

n{n + 1) 


[^70^^) - ^7^^^) + 2/?o] , 


4X/ 

n{n -|- 1) 



Tn ref. these transformations are presented in Bjorken x space to all orders in as- 


( 16 ) 
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Note also that 


^Clfg{n)jET — , A'y^'g{n)jET — ^'yNsi'^)MS ■ 

In (^) a superscript ”0” is used for the corresponding anomalous dimensions in the 
LO approximation. 

Note that the transformation of the coefficient functions and anomalous dimensions 
from the MS to the AB scheme is given by eqs. (p!5| - |^), in which the factor 2/?7,(n + l) 
should be replaced by 1/n. 


In each scheme the parton densities at Qo = 1 GeV^, as in 
in the form: 


13], were parametrized 


xAu^{x,Ql) = r]uAuX'^'^xUy{x,Ql) , 
xAdy{x,Ql) = r]dAdx'"’‘xdy{x,Ql) , 
xASea{x,Ql) = risAsx°'^xSea{x,Ql) , 
xAG{x, Ql) = rigAgX°‘<’xG{x, Ql) 


(18) 


where on R.H.S. of ( p!8D we have used the recent MRST unpolarized densities [^. The 
normalization factors Af are determined in such a way as to ensure that the first 
moments of the polarized densities are given by r//. 

The first moments of the valence quark densities rju and rjd are fixed by the octet 
nucleon and hyperon [3 decay constants ||^ 


gA = F + D = 1.2573 ± 0.0028, ag = 3F - D = 0.579 ± 0.025 . 
and in the case of SU(3) flavour symmetry of the sea (A-u = Ad = As at Ql) 

r]u = 0.918 , Tjd = —0.339 . 

The rest of the parameters in (^) 

dSi ^Vg 5 ! 


(19) 


( 20 ) 


( 21 ) 


have to be determined from the best fit to the A^{x, Q^) data. 

To calculate the virtual photon-nucleon asymmetry A^{x,Q‘^) in NLO QCD and 
then fit to the data we follow the procedure described in detail in our previous papers 


13], where the connection between measured quantities, Wilson coefficients and par- 


ton densities is given. 


5 


















The numerical results of the fits in the JET, AB and MS schemes to the present 
experimental data on |I|, ||, ||, |2^, ^ are listed in Table 1.^ The 

data used (118 experimental points) cover the following kinematic region: 

0.004 < X < 0.75, 1 < < 72 GeV"^ . (22) 


In this paper we present the results of the fit to the data averaged over 
at each x. (In our previous work |^) we also analyzed data in (x, Q^) bins. Our 
conclusion in the present paper also hold for this type of fit.) The total (statistical and 
systematic) errors are taken into account. The results presented in Table 1 correspond 
to an SU(3) symmetric sea. Note that in this case the first moment of the strange sea 
quarks r]s = As(1,Qq) = rjs/G . As was shown in our previous work |^, the flavour 
decomposition of the sea does not affect the quality of the fit and the results on the 
polarized parton densities AE, As and AG. 

I , the value of ttg is not well determined by the 


As in the case of the MS scheme 
fits in the JET and AB schemes to the existing data , i.e. /DOF practically does 

not change when Qg varies in the range:0 < < 1. In Table 1 we present the results 

of the fits corresponding to = 0.6. 


Table 1. Results of the NLO QCD fits in the JET, AB and MS schemes to 
the world data (Qq = 1 GeV'^). The errors shown are total (statistical 
and systematic). Qg = 0.6 (fixed). 


Scheme 

JET 

AB 

MS 

DOF 

118 - 5 

118 - 5 

118 - 5 


86.39 

86.35 

86.11 

XVDOF 

0.764 

0.764 

0.762 


0.267 ± 0.035 

0.267 ± 0.036 

0.255 ± 0.028 


0.124 ± 0.123 

0.124 ± 0.125 

0.148 ± 0.113 

as 

1.469 ± 0.460 

1.558 ± 0.583 

0.817 ± 0.223 

Vs 

- 0.027 ± 0.004 

- 0.022 ± 0.005 

- 0.049 ± 0.005 

Vg 

0.50 ± 0.12 

0.56 ± 0.14 

0.82 ± 0.32 

AE(1) 

0.416 ± 0.036 

0.444 ± 0.040 

0.287 ± 0.041 

ao(l GeV^) 

0.30 ± 0.04 

0.31 ± 0.05 

0.29 ± 0.04 


t After the completion of this work new data on Af and gf have been reported by the HERMES 
Collaboration |3^. 
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It is seen from the Table 1 that the values of /DOF coincide almost exactly 
in the different factorization schemes, which is a good indication of the stability of the 
analysis. The NLO QCD predictions are in a very good agreement with the presently 
available data on ^ as is illustrated in the JET scheme £t in Fig. 1. We would like 
to draw special attention to the excellent fit to the very accurate E154 neutron data 

= 1.1 for 11 experimental data points). 

The extracted valence and gluon polarized densities at Qg = 1 GeV^ for the different 
schemes are shown in Fig. 2 (solid, dotted and dashed curves correspond to JET, AB 
and MS scheme, respectively). Note that the valence densities xAuy and xAdy in the 
JET and AB schemes are almost identical so the dotted curves corresponding to xAuy 
and xAdy are not shown in Fig. 2. The difference between xAuy and xAdy in the MS 
and JET scheme is negligible. The results of the fit for the polarized valence densities 
are in an excellent agreement with what follows from the theory, namely, that they 
should be the same in the factorization schemes under consideration. 

The results on the polarized gluon densities determined by the fit in the JET, AB 
and MS schemes are also consistent. The values of their first moments rjg coincide 
within errors (see Table 1). However, as a consequence of the uncertainty in determin¬ 
ing the gluon density from the present data, the central values of r]g and therefore, the 
gluon densities themselves, differ somewhat in the various schemes. 

In Table 1 we also present our results in the different schemes for the hrst moments 
AE(1) of the polarized singlet quark density as well as for the axial charge oq. The 
obtained singlet densities AS(x, Ql) are shown in Fig 3a. It is seen from the table that 
the first moments AE(1) in the JET and AB schemes are in a very good agreement. 
(We recall that according to the definition of the JET and AB schemes AE(1) should 
be the same in the both schemes.) The corresponding densities A'E{x,Q'^), however, 
are slightly different (see Fig. 3a) because their higher moments are not equal. 

Our result for using the values of AE(l);^ and {rig)jjg from 

Table 1 and the NLO transformation rules (see eqs. (|^) and (|T^) for n = 1) 

~ 0.476 ± 0.084 (23) 

coincides within errors with the values of AH{1) jet and AE(1)^^ 

AE(1) JET = 0.416 ±0.036, AE(l)Aij = 0.444 ± 0.040 (24) 

determined directly by the £t to the data in the JET and AB schemes as presented 
in Table 1. The singlet densities AI]{x, Ql) jet (solid curve) and AT,{x,Ql)j^^j^rp 
(dashed curve) are shown in Fig. 3b. 
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Finally, we would like to draw attention to the excellent agreement between the 
values of the axial charge ao(l GeV'^') determined in the different schemes (see Table 
1), which illustrates impressively how our analysis respects the scheme-independence 
of physical quantities. 

In conclusion, we have performed a next-to leading order QCD analysis of the 
world data on inclusive polarized deep inelastic lepton-nucleon scattering in the JET, 
AB and MS schemes. The QCD predictions have been confronted with the data on 
the virtual photon-nucleon asymmetry , rather than with the polarized 

structure function , in order to minimize higher twist effects. Using the 

simple parametrization (^) (with only 5 free parameters) for the input polarized parton 
densities it was demonstrated that the polarized DIS data are in an excellent agreement 
with the pQCD predictions for A^{x, Q'^) in all the factorization schemes considered 

Moreover, we have demonstrated the consistency between the results of the analysis 
in each scheme and the NLO transformation equations relating them. Such consistency 
gives us conhdence that the NLO analysis is reliable. 
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Figure Captions 


Fig. 1 Comparison of our NLO results in the JET scheme for [x, Q^) with the ex¬ 
perimental data at the measured x and values. Errors bars represent the total error. 

Fig. 2 Next-to-leading order input polarized valence and gluon distributions at 
Q2 = 1 GeV^ in different factorization schemes. Solid, dotted and dashed curves 
correspond to the JET, AB and MS scheme, respectively. 

Fig. 3 (a) Next-to-leading order input polarized singlet distributions AE(x) at = 
1 GeV^ in different factorization schemes. Solid, dotted and dashed curves correspond 
to the JET, AB and MS scheme, respectively, (b) Comparison between the singlet 
density AT,{x)jet obtained from the fit (solid curve) and AT,{x)jjg^j^j. determined 
by eq. (H) (dashed curve), at = 1 GeV"^ . 
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